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ABSTRACT 

Using current observational data and simple dynamical modeling, we demonstrate that uj Cen is not 
special among the Galactic globular clusters in its ability to produce and retain the heavy elements 
dispersed in the AGB phase of stellar evolution. The multiple stellar populations observed in uj Cen 
cannot be explained if it had formed as an isolated star cluster. The formation within a progenitor 
galaxy of the Milky Way is more likely, although the unique properties of uj Cen still remain a mystery. 

Subject headings: globular clusters: individual (NGC 5139), Galaxy: formation 



1. Omega Centauri, the most massive globular 
cluster in the galaxy 

Omega Centauri is unique among globular clusters in 
the Milky Way in that its member stars exhibit a wide 
spread in metallicity (Freeman & Rodgers 1975; Freeman 
& Norris 1981; Mallia & Pagel 1981). This has led to 
the general belief that uj Cen was chemically self-enriched, 
while the other clusters were not (Suntzeff & Kraft 1996; 
see Meylan & Heggie 1997 for a review). In this context, 
it is often noted that uj Cen is the most massive globu- 
lar cluster in the Milky Way and thus possibly the most 
capable of retaining its own stellar ejecta (e.g., Ikuta & 
Arimoto 2000). The globular cluster Gl in the halo of the 
Andromeda galaxy has similar properties; it appears to be 
self-enriched and is also one of the most massive clusters 
in its host galaxy (Meylan et al. 2001). 

Omega Centauri has several other intriguing and po- 
tentially related properties. The metal-rich stellar pop- 
ulation is more centrally concentrated and has a lower 
velocity dispersion than the metal-poor population (Nor- 
ris et al. 1997). Moreover, the metal-rich population has 
no detectable rotation, while the metal-poor population 
rotates relatively rapidly (~ 7 km s _1 ), consistent with 
the apparent flattening of the cluster, which is among the 
highest of all globular clusters in the Galaxy (Frenk & 
Fall 1982). The metal-rich stellar population also appears 
to be younger than the metal-poor population, consistent 
with the idea of self-enrichment over a period At = 3 — 5 
Gyr (e.g., Hughes & Wallerstein 2000; Hilker & Richtler 
2000; Lee et al. 1999). This is supported by the large 
enhancement in s-process elements in the metal-rich pop- 
ulation relative to that in globular clusters with similar 
metallicities (Norris & Da Costa 1995; Vanture, Waller- 
stein & Brown 1994). These properties have led to the 
suggestion that uj Cen managed to retain the ejecta from 
the first generation of stars in the asymptotic giant branch 
(AGB) phase, which then formed a centrally concentrated, 
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metal-enriched, younger generation of stars (e.g. Norris et 
al. 1996; Smith et al. 2000). 

In this paper we demonstrate that uj Cen is not special 
in the ability to retain its own stellar ejecta. The escape 
speed in uj Cen is not the highest among the Galactic glob- 
ular clusters. We find that more than half of the clusters 
could retain at least 20% of their AGB winds. It is there- 
fore surprising that uj Cen is the only globular cluster in 
the Milky Way to show signs of self-enrichment. 3 We con- 
clude that this cluster could not have evolved in isolation 
on its present orbit in the Galaxy and that it must have 
been substantially different in the past. We discuss other 
models suggested previously, including the idea that uj Cen 
is the nucleus of a dwarf galaxy that lost its outer envelop 
by tidal stripping (Zinnecker et al. 1988; Freeman 1993). 

2. THE ESCAPE VELOCITY FROM THE CLUSTER 

To compute the escape velocities of Galactic globular 
clusters, we use the photometric data from the catalog 
by Harris (1996). Out of 147 clusters in the catalog, we 
choose 141 for which the data are available on the distance, 
integrated magnitude, core radius, and concentration pa- 
rameter. All clusters are fit by single-mass isotropic King 
models with a constant mass-to-light ratio, MjLy = 3 in 
solar units. 

First, we check that the photometric models are consis- 
tent with the directly measured central velocity dispersions 
for a sample of 56 clusters from Pryor & Meylan (1993), 
their Table 2. The most reliable estimates of the disper- 
sion are based on individual stellar radial velocities, but 
these estimates depend on the position of stars relative to 
the cluster center. The core velocity dispersion of uj Cen 
is 17 km s _1 (Merritt, Meylan & Mayor 1997), with the 
indication that it rises in the center up to 25 km s" 1 (van 
Leeuwen et al. 2000). 

We thus avoid a direct comparison and show instead in 
Figure 1 two separate histograms of the central velocity 
dispersions as inferred from the photometric models and 
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as measured directly. The distributions are qualitatively 
similar. Also, the photometric model of uj Cen predicts a 
central velocity dispersion of about 17 km s _1 , the same 
as the observed value. 



sition: 

fw = J f( E ) drdv ^ J dn u h(v csc (r) - |v + uj), 

(2) 

where Cl u is a solid angle of the wind velocity u w , and h(x) 
is a unit step function. 
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Fig. 1. — Distribution of the central line-of-sight veloc- 
ity dispersion of the Galactic globular clusters, as inferred from 
the photometric models (left panel) and observed spectroscopi- 
cally (right panel). The arrows point to the velocity dispersion 
of lo Cen. 

We note that in the past 20 years there has been only 
one study (Dubath et al. 1997) providing homogeneous ve- 
locity dispersion measurements for a significant number of 
clusters. Now, with the availability of large-field survey 
spectrographs such data should be easier to obtain and 
would provide a very useful constraint on the mass models 
of globular clusters. 

Having calibrated the photometric models, we can use 
them to determine the escape velocity. This is calcu- 
lated using the reduced gravitational potential, w csc = 
(2$ t — 2$) 1 / 2 , in the center (giving the maximum escape 
velocity) and at the half-mass radius (providing a thresh- 
old for most stars). Here, the tidal potential <T> t accounts 
for truncation at the tidal radius. 

Figures 2 and 3 demonstrate that u Cen is not special 
in its escape velocity. There arc 11 clusters with higher 
v OS c,o an d 8 clusters with higher 

^esc,h- 

We also find that the ratio of the central escape velocity 
to the central velocity dispersion can be fit by a simple 
relation 



: ,o/<7o = 3.7 + 0.9 (c- 1.4), 



(1) 



where c is the concentration of the King models. The fit 
is accurate to 3% in the range c = 0.4 — 2.8. 

We now calculate the fraction of AGB winds retained 
by a cluster as a function of the wind velocity u w . We 
integrate the distribution function of stars, f(E), over the 
phase space volume within which the total velocity of the 
wind particles is below the escape speed at the star's po- 
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Fig. 2. — Central escape velocity for 141 Galactic glob- 
ular clusters vs their mass M. Omega Centauri is denoted by 
an asterisk. The escape velocity is calculated from the photo- 
metric single-mass King models with M/L v = 3. 
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Fig. 3. — Escape velocity at the half-mass radius for the 
Galactic globular clusters vs their mass M. Omega Centauri is 
denoted by an asterisk. 

Figure 4 shows the result for clusters of various concen- 
tration parameters, from c = 0.6 to 2.4. The fraction f w 
declines gradually and even the clusters with u w ~ 3co 
can retain 10 to 20% of their winds. Since according to 
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Loup ct al. (1993) the characteristic terminal velocity of 
AGB winds is u w w 15 km s _1 (although a velocity as 
high as 90 km s _1 has been measured for a field red giant 
by Dupree, Sasselov & Lester 1992), more than half of the 
Galactic clusters should be able to retain a sizable fraction 
of their winds. 
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Fig. 4. — Fraction of AGB winds retained by clusters 
of concentration c (marked by numbers) as a function of the 
wind velocity in units of the central velocity dispersion. The 
top solid line shows the model with the highest binding en- 
ergy (c = 1.4), the dotted (c = 0.6) and dashed (c = 2.4) 
curves give the lower bounds. All models with intermediate 
concentrations fall between the plotted lines. For all con- 
centrations, the fraction can be fit reasonably accurately by 
f w (x = u w /a ) = 1/[1 + (x/2f + (x/4) 4 + (x/3) 8 }. 



3. THE ORBIT AND THE POSSIBILITY OF 
SELF-ENRICHMENT 

The question now is: could lo Cen (and other clusters) 
retain the wind material throughout their evolution in the 
Galaxy? The oldest stars in the Galactic disk are dated to 
be at least 9.5 ± 0.8 Gyr old (Oswalt et al. 1996; Leggett 
et al. 1998). The age of the Universe is determined fortu- 
itously well from the coincidental CMB correlation (Knox 
et al. 2001) to be 14.0 ±0.5 Gyr. Thus, it is quite possible 
that the dense gaseous disk was already assembled when 
the Universe was about 4 Gyr old, equal to the age spread 
of stars in lo Cen. If this was the case, ram pressure of the 
disk would strip any gas left in lo Cen, as we show below. 

The current Galactocentric velocity of lo Cen from the 
recently measured proper motion (Dinescu et al. 1999) is 
relatively low: Vr = —31 km s _1 , Vq = —65 km s _1 , 
V z — 4 km s _1 , with the uncertainty ±10 km s _1 . The 
speed of only 70 km s _1 makes the cluster strongly bound 
and thus, its present position at 6.3 kpc from the Galac- 
tic center is likely to be close to the orbital apocenter. 
Indeed, calculations of Dinescu et al. (1999) give the peri- 
center distance R p = 1.2 kpc and the apocenter distance 
R a = 6.4 kpc. Therefore lo Cen remains close enough to 
the Galactic center to pass through the disk twice in each 



orbital period, which is relatively short, P or b = 1.2 x 10 8 

yr- 

We calculate the amount of gas that might have been 
present in lo Cen as follows. The total mass-loss rate of 
AGB stars depends on the slope of the IMF between 1 
and 8 M©, dN/dm oc m~ 7 , on the mass of the white 
dwarf remnant, m w d, and on the main sequence lifetimes 
of stars, t ms : 

dM w dN ( dtmsY 1 

— = d^ {m - mwd) {~^r) ' (3) 

where m = m ms (t) is the inversion of the function t ms (m). 
We consider two cases, 7 = 2.2 (following Kroupa, Tout & 
Gilmorc 1993) and 7 = 3.3 (Scalo 1986). Using the code 
of Hurley, Pols & Tout (2000), we calculate t ms {m) for the 
oldest stars in lo Cen, with the metallicity [Fe/H] = —1.7, 
and find that in the mass range of interest it is well fit 
by t ms w 5 x 10 9 (m/M Q )~ 2 - 4 yr. The masses of stel- 
lar remnants have been determined observationally in the 
nearby star clusters (Weidemann 2000) , and the result can 
be accurately fit by m w d = 0.444 M© + 0.084 m ms . Com- 
bining all these expressions, we find the mass- loss rate 
at the age of 4 Gyr: dM w /dt « 10" 2 M cX Gyr" 1 . Here 
M c i ps 3 x 10 6 Mq is the present mass of the cluster. (In- 
terestingly, the mass-loss rate is almost the same for both 
values 7 = 2.2 and 7 = 3.3.) Thus, in half the orbital 
period the cluster would accumulate M g = 2 x 10 3 M Q of 
gas. 

We integrate the orbit of lo Cen backward in time using 
an accurate Galactic potential that fits most observational 
constraints of the Milky Way (model A\ of Klypin, Zhao & 
Somervillc 2002). We model the evolution of the Galaxy 
starting with a spherical halo 14 Gyr ago, compress the 
baryons into the disk and the bulge over an exponential 
timescale of 3 Gyr and then turn them into stars over the 
next 3 Gyr. The bulge is gradually converted into the 
present flattened shape. 

We find the locations where the cluster crosses the plane 
of the disk and then compare the ram pressure of the 
disk gas with the restoring gravitational force per unit 
area at the half-mass radius of lo Cen. In many passages 
PdiskWroi > GM c iM g /(4nRf l ) and therefore the AGB ejecta 
should be stripped from the cluster. Also, the current mass 
of lo Cen is too small for dynamical friction to affect the 
orbit. 

To summarize, we have shown that if (1) lo Cen had 
formed in isolation with its present mass, and (2) the 
Galactic disk was in place prior to the formation of the 
metal-rich stars, then the cluster would encounter the disk 
gas in its orbital motion and lose the gas accumulated from 
the enriched AGB winds. The same argument applies to 
all globular clusters within 10 kpc from the Galactic cen- 
ter. They could not have formed stars from the material 
enriched by s-process elements. 

Alternatively, if all the stars in the oldest globular clus- 
ters form prior to the assembly of the Galactic disk, then 
lo Cen must be unrealistically special. Why did it retain 
its enriched gas and formed new populations of stars over 
4 Gyr whereas a few dozen other clusters, capable of doing 
so, did not. As we have shown, even though lo Cen is the 
most massive cluster its potential well is similar to many 
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others. So we are left again with the puzzle of why u> Cen 
is different from all other Galactic globular clusters. 

We conclude that it is highly unlikely that u Cen formed 
and evolved in isolation and enriched itself with heavy el- 
ements of the AGB winds. 

4. IN CLOSING: TOWARDS A CORRECT STORY OF UJ CEN 

In light of the unusual properties of u> Cen, it is becom- 
ing a popular scenario that the cluster may be the core 
of a disrupted dwarf galaxy (e.g., Freeman 1993). The 
attraction is that a continuous infall of gas to the center 
of the dwarf galaxy may lead to a variable star formation 
history. Kormendy (1985) has shown that dwarf galaxies 
split on the plot of the central surface brightness vs core 
radius into dense dE galaxies (like M32) with the central 
velocity dispersion a in excess of 50 km s -1 , and diffuse 
dSph galaxies (like Draco) with ao < 20 km s _1 . Tidal 
stripping of the outer envelope would decrease the total 
luminosity of the dwarf galaxy but would not significantly 
change its core properties, /xo and R c . Thus, the cores of 
dE would have too high a velocity dispersion and the cores 
of dSph would have too low density to count as globular 
clusters. 

The counter-examples to this rule are nucleated dSph 
galaxies, sometimes classified as dE,N (Peterson & Cald- 
well 1993). NGC 205 (Meylan 2001) has a core velocity 
dispersion of 30 km s _1 , dropping to 15 km s _1 in the 
bright compact nucleus (Carter & Sadler 1990). The re- 
cently discovered Sagittarius galaxy apparently has a glob- 



ular cluster M54 positioned at its center (Ibata, Gilmore 
& Irwin 1994). Also, Drinkwater et al. (2000) detected 
five Ultra-Compact Dwarf galaxies in the Fornax cluster 
of galaxies, which are brighter than globular clusters but 
almost as compact. They could be the result of mergers of 
young star clusters, as suggested by Fellhauer & Kroupa 
(2002). It remains to be seen if the nuclei of these galaxies, 
when stripped of their envelopes, would have similar kine- 
matic and chemical properties to u> Cen. Also, detailed or- 
bital calculations are necessary to determine if the cluster 
can be brought in to its present location near the Galactic 
center (see Zhao 2001). 

The uniqueness of u> Cen among the Galactic globular 
clusters still poses unanswered questions. In the hierarchi- 
cal formation scenario we expect a few dozen progenitor 
galaxies to make up the Milky Way. If u> Cen did in fact 
form at the center of one of them, why did not other al- 
most as massive clusters form in the other progenitors and 
have similar chemical enrichment? Were they entirely dis- 
rupted or was the history of ut Cen completely different? 
We intend to address these issues in forthcoming work. 

We thank Chris Tout for discussions. OG is supported 
by the STScI Institute Fellowship, and JEP acknowledges 
the continuing support by the STScI visitors program. 
The electronic tables of the escape velocities, velocity 
dispersions, and main-sequence lifetimes are available at 
http : //www . stsci . edu/^ognedin/gc/. 
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